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The  fuel  cell  technology  is  a  key  element  for  the  hydrogen  energy  economy  and  therefore  crucial  for 
sustainable  development.  High  temperature  proton  exchange  membrane  (HT-PEM)  fuel  cells  (FC)  can  be 
operated  with  reformate  gas  and  thus  represent  an  important  bridging  technology  for  the  energy 
transition  to  a  renewable  energy  based  system.  HT-PEM  FCs  based  on  phosphoric  acid  (PA)  are  still 
subject  to  intense  research,  investigating  the  electrolyte  behaviour.  By  enhancing  state  of  the  art  2D  FEM 
simulations  of  FCs  with  the  vapour  liquid  equilibrium  of  water-phosphoric  acid  and  evaporation  ki¬ 
netics,  a  model  was  created  in  which  the  local  concentration  of  PA  can  be  calculated.  Knowledge  of  the 
concentration  field  yields  the  basis  for  calculating  the  locally  varying  ionic  conductivity  and  other 
physical  properties.  By  describing  the  volume  expansion  behaviour  of  PA  it  was  possible  to  predict  the 
catalyst  particle  deactivation  due  to  the  swelling  of  PA.  The  in  situ  concentration  predicted  by  the 
simulation  ranges  from  96  to  111  wt%.  The  model  was  validated  using  measured  data  of  a  single  cell 
design  for  different  temperatures  and  pressures.  By  varying  the  PA  content  flooding  of  the  simulated  fuel 
cell  could  be  observed  and  was  linked  to  humidification  effects. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  efficient  energy  converters  that  convert  the 
chemical  energy  of  a  fuel  into  electrical  energy.  In  addition  to  high 
conversion  efficiency  compared  to  conventional  technology  the  key 
advantage  of  fuel  cells  is  the  possibility  of  a  renewable  fuel.  The 
usage  of  hydrogen  from  renewable  sources  yields  zero  emissions, 
local  as  well  as  global. 

The  market  shows  rapidly  growing  interest  in  the  transport 
sector.  For  stationary  applications  the  first  product  lines  are  already 
available:  Fuel  cells  equipped  with  internal  gas  reformers 
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generating  both  power  and  heat.  The  mentioned  application  as  a 
micro  combined  heat  and  power  system  is  of  great  importance  for 
Germany  [1],  The  operation  of  the  micro-CFIP  with  substitute  nat¬ 
ural  gas  from  renewable  sources  is  one  of  the  first  major  milestones 
in  the  energy  transition.  The  usage  of  renewable  synthetic  gas  re¬ 
quires  gas-to-power  technologies  and  places  high  demands  on  the 
flexibility  of  the  power  converter,  especially  with  respect  to  the 
fuel.  HT  PEM  FCs  are,  because  of  their  high  CO  tolerance  [2],  suitable 
for  use  with  Fl2-rich  synthetic  natural  gas.  Furthermore,  they  also 
achieve  excellent  efficiency. 

The  critical  component  of  a  FC  is  the  membrane  electrode  as¬ 
sembly.  Recent  research  focuses  on  phosphoric  acid  doped  poly¬ 
benzimidazole  (PBI)  membranes  for  HT  PEM  FC  applications  [3,4], 
These  have  the  advantage  of  being  ionically  conductive  at  low 
relative  humidity  and  therefore  auxiliary  equipment  regarding 


T.J.  Kazdal  et  al.  /  Journal  of  Power  Sources  249  (2014)  446—456 


447 


Nomenclature 

Mi 

molar  mass 

Na 

Avogadro’s  number 

kB 

Boltzmann  constant 

P 

pressure 

a 

transfer  coefficient 

Pvp 

vapour  pressure 

AH° 

standard  reaction  enthalpy 

Q 

heat  source/sink 

8 

diffusion  length 

R 

gas  constant 

Pact 

activation  polarization 

Ri 

species  source/sink 

T 

pressure  coefficient 

Rcat 

electrode  roughness 

Kbr 

permeability 

ref 

reference 

A; 

stoichiometric  coefficient 

T 

temperature 

®ik 

binary  diffusion  coefficient 

u 

velocity 

A 

dynamic  viscosity 

vm 

molar  volume 

w 

mass  fraction 

X 

doping  level 

0 

volume  fraction 

Xi 

mole  fraction 

P 

density 

A 

area 

a 

conductivity 

carb 

carbon  black 

X 

tortuosity 

CL 

catalytic  layer 

©active 

ECSA  per  volume  catalyst 

ECSA 

electrochemical  active  surface  area 

©incorp 

incorporated  platinum  area 

GDL 

gas  diffusion  layer 

e 

porosity 

HOR 

hydrogen  oxidation  reaction 

^electrode 

,  electrode  thickness 

J 

current  density 

q 

concentration 

m 

mass 

cp 

heat  capacity 

MEA 

membrane  electrode  assembly 

Dik 

diffusion  coefficient 

No. 

number  of 

dpt 

diameter  platinum  particle 

ORR 

oxygen  reduction  reaction 

E 

electrode  potential 

PA 

phosphoric  acid 

Ea 

activation  energy 

PAM 

phosphoric  acid  mixture 

F 

Faraday  constant 

PBI 

polybenzimidazole 

Hi 

Henry  coefficient 

PC 

phase  change 

io 

exchange  current  density 

STP 

standard  conditions  for  temperature  and  pressure 

k 

thermal  conductivity 

TGA 

thermogravimetric  analysis 

feat 

catalyst  loading 

humidification  is  assumed  to  be  omittable  [5],  However,  many 
physical  properties  of  phosphoric  acid,  as  well  as  the  ionic  conduc¬ 
tivity,  are  a  function  of  the  water  content  in  the  acid  [6—8].  Unlike  for 
low  temperature  (LT)  PEM  FCs,  for  the  high  temperature  (HT)  PEM 
FCs  only  very  few  studies  on  the  water  balance  were  performed. 

In  the  literature,  many  approaches  to  model  HT  [9,10]  PEM  FCs 
[11,12]  can  be  found.  Until  now,  the  distribution  of  phosphoric  acid 
was  not  subject  to  the  investigations.  Within  this  paper  the  2D 
model  of  a  HT  PEM  FC  with  a  PBI  membrane  electrode  assembly 
(MEA)  is  presented.  Special  emphasis  lies  on  the  water  balance  in 
the  membrane.  Understanding  the  balance  is  crucial  for  the  proper 
simulation  of  HT  PEM  FCs  because  the  assumption  of  constant, 
water  independent  physical  properties  is  not  necessarily  valid. 
With  a  model  for  the  vapour  liquid  equilibrium  the  water 
household  of  the  FC  was  calculated.  On  the  basis  of  the  obtained 
distribution  of  water,  the  physical  parameters  of  the  electrolyte 
were  determined.  Moreover,  the  sorption  behaviour  of  PA  was 
investigated  and  its  volume  expansion  and  thus  the  pore  filling 
degree  was  implemented  into  the  model.  Since  a  flooded  pore  is 
assumed  to  be  electrochemically  deactivated,  the  reaction  kinetics 
were,  in  consequence,  coupled  with  the  local  concentration  of 
water. 

2.  Model 

2.1.  Model  description  and  assumptions 

In  this  work  a  2D  model  of  a  high  temperature  PEM  fuel  cell  is 
presented.  The  calculation  domain  consists  of  gas  channels,  gas 


diffusion  layers,  catalytic  layers  and  a  membrane,  as  seen  in  Fig.  1. 
The  dimensions  are  summarized  in  Table  1.  The  following  as¬ 
sumptions  were  made: 

2.1.1.  General 

•  Anisothermal 

•  Steady  state 

•  ORR  dominates  activation  overpotentials 

•  No  gas  permeation  through  the  membrane  — >  no  mixed 
potential 

•  Operation  on  hydrogen  and  air 

•  Ideal  gas  mixtures  in  gas  phase 


lation  domains. 
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Table  1 

Model  dimensions. 


H2  channel 
Anode  GDL 

Cathode 
Cathode  GDL 
Air  channel 
Channel  depth 
Channel  length 


2.3.2.  Materials 


enters  and  leaves  the  membrane.  R.  He  [17]  studied  the  influ¬ 
ence  of  the  doped  acid  and  absorbed  water  on  the  volume 
swelling  for  acid  doped  PB1  membranes  and  found  the  effect  of 
swelling  by  water  to  be  negligible 

•  Ionic  conductivity  and  diffusion  coefficient  depend  on  acid 
concentration,  temperature  and  morphology  (e,  x) 

•  Water  gradient  through  the  membrane  is  of  higher  magnitude 
than  along  the  channel 

In  Fig.  2  the  magnified  model  geometry  is  shown  and  the 
naming  convention  for  the  domains  and  the  boundaries  is  pre¬ 
sented.  The  governing  equations  and  their  computational  domains 
are  summarized  in  Fig.  3. 


•  Macro  homogeneous  and  isotropic 

•  Pore  size  distribution  in  the  catalytic  layer  has  two  maxima 

•  Effect  of  porosity  on  transport  parameters  is  described  by  the 
Bruggemann  correlation 

2.3.3.  Phosphoric  acid 

It  is  assumed  that  phosphoric  acid  polymerises  as  described 
by  R.F.  Jameson  [13],  In  his  work  he  presents  the  equilibrium 
concentrations  of  the  polyphosphoric  acids  at  15  °C.  To  account 
for  the  temperature  dependence  of  the  reaction  equilibrium  we 
used  data  from  R.A.  Munsen  [14,15],  He  studied  the  kinetics  of 
the  ortho-pyro  interconversion  and  found  the  equilibrium  con¬ 
stant  to  be  of  the  magnitude  of  E-03  in  the  range  of  40-100  °C. 
Therefore  it  is  reasonable  to  neglect  the  temperature  depen¬ 
dence  and  to  use  Jameson’s  data  for  the  estimation  of  the 
composition  distribution  while  discussing  the  concentration  of 
PA. 

Furthermore  for  the  physical  properties  of  phosphoric  acid  the 
effects  of  this  polymerization  are  accounted  for  by  the  imple¬ 
mentation  of  the  concept  of  apparent  concentrations:  A  mixture  of 
90  wt%'  phosphoric  acid  heated  to  200  °C  in  a  closed  system  has  a 
metrological  determinable  ionic  conductivity.  Its  value  is  related  to 
the  temperature  and  the  apparent  concentration  of  90  wt%  PA, 
despite  the  fact  that  polycondensation  might  have  occurred.  This 
concept  was  applied  to  all  empirically  determined  physical 
properties. 

•  no  polycondensation  thermodynamics  and  kinetics  of  phos¬ 
phoric  acid 

-  Only  one  species  ->  concept  of  apparent  concentration 

—  Physical  properties  (D,<j,pvp)  still  account  for  poly¬ 
condensation  and  the  actual,  but  unknown  composition 

•  Diffusion  coefficient  of  PA  depends  on  temperature  and 
concentration 

•  Vapour  pressure  of  species  PA  is  neglected.  At  221  °C  E.H. 
Brown  [16]  observed  a  mass  fraction  of  0.9  wt%  in  the  vapour 
phase. 

•  Acid  expansion  and  contraction  is  completely  reversible  during 
one  simulation  run 

•  Capillarity  does  not  affect  evaporation 

2.3.4.  PBI-membrane 

•  Contains  free  acid  and  can  be  considered  as  a  porous  medium 

•  No  thickness  change  over  lifetime 

•  No  volume  swelling  during  a  simulation  run.  Assuming  that  the 
amount  of  PA  is  constant  during  a  simulation  run,  only  water 


2.2.  Governing  equations 

2.2.3.  Conservation  of  mass  and  momentum 

The  velocity  field  is  obtained  by  solving  the  equation  of  conti¬ 
nuity  in  combination  with  the  Navier- Stokes  and  the  Brinkman 
equations  for  flow  in  porous  media.  It  is  important,  to  mention  the 
porosity  dependence  of  Brinkman’s  equation,  as  seen  in  EQ:1, 
because  Chapter  2.3  will  introduce  the  porosity  in  the  catalytic 
layer  as  a  dependent  variable. 

c(U'VQ)  =V[“P£+c(VU  +  (VU)T)“lc(VU)£]  “(^r  +  rh)U 

(1) 

There  are  two  convective  mass  flows  to  consider,  one  in  each 
channel.  The  gas  diffusion  layers  (GDL,  [2]  &  [6])  and  the  catalytic 
layers  (CL,  [3]  &  [5])  are  modelled  as  porous  media.  The  energy 
balance  and  the  conservation  of  species  are  coupled  with  the 
conservation  of  mass,  to  account  for  the  change  of  density  and 
temperature.  Table  2  shows  parameters  for  calculating  the  con¬ 
servation  of  mass  and  momentum. 

2.2.3.3.  Boundary  conditions,  uhor  and  uorr  are  flow  velocities.  For 
stoichiometric  operation  the  velocities  depend  on  the  current 
generated  and  the  cell  geometry  according  to  EQ:2: 

Uhqr  = _ 1 _ Jhor  A>1)  •f?-r-AAno  . 

NO. channels  ’^Channel  2F  pH0K 


2.2.I.2.  Source  terms.  At  both  electrodes  mass  is  produced  and 
consumed.  The  two  effects  considered  are  the  reaction  conversion 


100%  ortho-PA  @  15  °C  [13], 


Fig.  2.  Model  geometry  and  naming  i 
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Navier-Stokes  — 
Maxwell-Stefan  — 
Maxwell-Stefan  — 
Ohm's  Law- 
Ohm's  Law  ~ 
Heat  Transfer  — 
Fick's  Law  — 


, 

: 

■ 

J 

t 

Anode  side 

Fig.  3.  Governing  equations  and  th 


Cathode  side 
r  calculation  domains. 


and  the  phase  change.  Reaction  conversion  is  directly  related  to  the 
electrical  flux,  as  seen  in  EQ:3  and  EQ:4.  The  change  of  mass  due  to 
phase  change  will  be  introduced  in  Section  2.3. 


(MH20  M0l\ 


2.2.2.  Conservation  of  species 

The  Stefan-Maxwell  equation  (EQ.:6)  is  used  to  model  the 
conservation  of  species.  By  solving  the  equations  for  the  conser¬ 
vation  of  mass  and  momentum  the  convective  velocity  u  is  pro¬ 
vided.  The  density  within  gaseous  phases  is  calculated  by  the  ideal 
gas  law.  For  the  calculation  of  the  Stefan-Maxwell  diffusion  co¬ 
efficients  for  a  three  component  system  EQ:7  [18]  can  be  used.  The 
binary  diffusion  coefficients  are  calculated  by  EQ:8  [18]  and  EQ:9 
[19]  respective  their  validity2  and  corrected  by  the  Bruggeman 
equation  (EQ:11)  [20]  regarding  the  porosity  of  the  respective 
domain  material. 


Vjj  +  p(uV)w;  =  R, 

if  =  -  (pot  Dik^k) 

(Oi(uk+O)m)  |  Wk(Uj  +  Um)  (d2 
XfOfen  XkVim  Xm®ifc 

lk~  *i  {  xk  |  xm 

®ik®im  ®ik®km  ®im®km 

0.00143  T175  (  1 


(5) 


Pc,iPc,k)  {TC,{TC 


2.2.2.I.  Sources.  The  source  terms  are  equal  to  those  presented  in 
Section  2.2.1.  For  the  conservation  of  species  they  are  mapped  to  the 
reacting  species.  E.g.  the  mass  loss  in  the  anode  is  a  loss  of  the 
species  hydrogen  (Table  3). 


binary  diffusion  coefficients  of  m 


Domain 

Important  parameter 

[2] 

e[2]  =  0.4 

[3] 

=  e(T,WpA) 

[5] 

=  £(T,6)pa) 

[6] 

e[6]  =  0.4 

Domain 

Initial  value 

[1] 

UHOR 

17] 

HI, [21, [3] 

Phor 

[5], [6], [7] 

PORR 

Boundary 

Condition 

14] 

Phor 

[B] 

[C] 

PI 

iihor.  laminar  inflow 

Source 

Term 

[3] 

UIhoR  +  rilpc 

[5] 

moRR  +  mpc 

2.2.3.  Conservation  of  charge 

2.2.3.I.  Electron  flux.  The  approach  for  describing  the  electrical 
flux  can  be  divided  into  two  regions:  anode  and  cathode.  Each 
region  has  an  electron  source  and  a  boundary  condition  where  the 
potential  is  constant,  creating  a  potential  field  over  the  domain. 
Combining  this  potential  field  with  the  electrical  conductivity 

leads  to  an  electron  flux.  As  seen  i 

n  ble  4  the  potential  field  is 

introduced  between  the  two  interior  boundaries  [a]  and  [b].  Cur¬ 
rent  sources  are  the  two  catalytic  layers  [3]  and  [5]  with  the 
hydrogen  oxidation  reaction  as  the  source  for  electrons  and  ions 

and  the  oxygen  reduction  reaction  a 

is  the  sink  for  both  species.  We 

use  the  Bruggeman  equation  [20] 

to  correct  the  electrical  con- 

ductivity  of  the  catalytic  and  gas  diffusion  layers  with  respect  to 
their  carbon  volume  fraction.  Electrons  are  produced  in  the  cata¬ 
lytic  layer  of  the  anode  [3]  and  with  the  lower  potential  at  [a]  we 
obtain  the  electrical  flux  for  the  anode  side.  At  the  cathode  the 
charge  travels  from  the  place  of  higher  potential  [b]  to  the  sink  at 
the  catalytic  layer  of  the  cathode  [5],  where  the  reaction  consumes 

them  (Tables  5). 

j  =  a E  +je 

(10) 

Parameter,  boundary  and  initial  values  regarding  the  conserva- 

tion  of  species. 

Domain 

[11,  [2],  [3] 

ojH2  =  1  -  ojH2q 

[5],  [6],  [7] 

Wq2  =  1  —  0JNl  -  Wh20 

Boundary 

Condition 

[4] 

k  =  o 

[B] 

k  =  o 

[C] 

«H2  =  1  -  <l>h20 

PI 

Wq2  =  1  —  0JNl  -  WH20 

Source 

Term 

[3] 

Bin,  =  Bihor 

[5] 

mo,  =Jorr(^) 

B1h]0  =  Jorr  +  mj^ 
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2.2.32.  Ion  flux.  Ion  flux  is  modelled  similar  to  electron  flux.  Two 
current  sources  ([3]  and  [5])  lead  to  a  potential  field  and  thus  a  flux 
of  ions.  The  potential  depends  on  the  ionic  conductivity  of  the 
electrolyte.  The  doped  PBI  membrane  is  considered  to  be  a  porous 
medium  filled  with  PAM.  Conductivity  is  modelled  as  a  heteroge¬ 
neous  mixture,  which  can  be  described  by  a  method  developed  by 
D.A.G.  Bruggeman  [20]: 

°rnix  =  4 '  f7cond.</,cond.  (H) 

where,  (rmjX  is  the  conductivity  of  the  mixture,  cond.  refers  to  the 
conductive  material  and  </>  is  the  volume  fraction  of  that  material 
in  the  mixture.  A  and  B  are  coefficients.  The  two  cases  to  discuss 
are  the  conductivity  within  the  membrane  and  the  catalytic 
layers: 

1.  In  the  membrane  the  electrolyte  is  the  conductive  phase.  Its 
volume  fraction  can  be  calculated  from  the  doping  level  of 
the  membrane,  as  seen  in  EQ:12.  Knowing  the  volume  frac¬ 
tion  and  the  conductivity  of  PAM  [6,8]  it  is  possible  to  fit 
EQ:11  on  literature  data  [21—24],  obtaining  the  coefficients 
presented  in  EQ.:13,  for  a  doping  level  dependent  membrane 
conductivity. 


M  (X-1.25)Vm,PA 

0PAM  ~  x-v+V. 


(12) 


ff!onic  =  °-42  •  <WT,  Wpa)  •  Cm1 


(13) 


2.  The  volume  of  the  catalytic  layers  ([3]  and  [5])  is  shared  by  the 
solid  electrode  material,  the  electrolyte  and  a  void  fraction.  In 
this  case  the  ion  conductivity  of  the  electrolyte  is  assumed  to 
follow  the  simplified  Bruggeman  equation,  EQ:  14: 


Table  6 

Initial  conditions  and  parameters  for  EQ:  16-19. 
Variable  Description 

a  Transfer  coefficient 

x,p  Partial  pressure  reactant  i 

pret  Reference  pressure 

y  Pressure  coefficient 

Ea  Activation  energy 

Tef  Reference  temperature 

Icat  Catalyst  loading 

dpt  Diameter  platinum  particle 

Beiectrode  Electrode  thickness 


Value  (Ano.|Cath.) 


0.5 1 1 

(ORR:  67  kj  mol-1) 
25"  C 

0.7 1 1  mg  cm-2 

50|60  nm 
0.995 


•fiS  =  ffPAMCT>WRA)-0PAM  (14) 

2.2.4.  Conservation  of  energy 


pCpUVT  =  V-(kVT)  +Q.  (15) 

The  energy  balance  is  modelled  thoroughly,  including  evapo¬ 
ration  and  condensation  enthalpies.  Due  to  the  conditions  that 
incoming  fluxes  and  walls  are  kept  at  a  constant  temperature  the 
maximum  temperature  variation  experienced  was  below  7  K  for 
the  operating  voltage  of  0.35  V  at  433  K.  It  is  assumed  that  the  effect 
of  a  non  uniform  temperature  distribution  will  become  more 
important  for  3D  models  of  the  entire  cell. 

2.2.5.  Butler—  Volmer  equation 

The  Butler-Volmer  equation  describes  how  the  electric  current 
changes  with  the  electrode  potential.  Most  of  the  source  and  sink 
terms  introduced  earlier  depend  on  the  current  density  as 
described  by  the  Butler-Volmer  equation: 

(.6, 

The  modelling  approach  for  the  exchange  current  density,  io,  is 
derived  from  Barbir  [25],  The  reference  exchange  current  density  is 
from  Weber  [26]  and  consistent  with  Huang  [27]  and  Cheddie  [10], 
who  used  experimental  data  from  Parthasarathy  et  al. 


Important  parameter 


Parameter,  boundary  and  initial  values  for  EQ:20. 


Important  parameter 


[3],  [4],  [5] 
Source 
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^  ~  PptdptBeTectrode  0  ®in“”p)  ‘  ' '  ’ ' 

Rcat  is  called  electrode  roughness  and  connects  the  micro  geo¬ 
metric  dimensions  of  the  catalyst  particle  with  the  macroscopic 
dimensions  of  the  catalytic  layer  [25,10],  The  total  surface  area  of 
platinum  is  6Lcat/PptdptBeiectrode  if  spherical  particles  of  diameter 
dpt  are  assumed.  Its  unit  is  area  platinum  per  geometric  volume 
catalytic  layer,  ©incorp  is  used  to  account  for  platinum  incorpo¬ 
ration  and  deactivation.  Its  value,  being  close  to  unity  and 
therefore  considerably  reducing  the  active  surface  area,  corre¬ 
sponds  to  the  idea  that  only  a  small  part  of  the  catalytic  layer  is 
active  [28],  The  change  of  the  electrochemical  active  surface  area, 
©active,  was  introduced  as  a  relative  quantity  ranging  from  zero  to 
unity.  It  modifies  the  exchange  current  density  as  phosphoric 
acid  penetrates  the  catalytic  layer.  It  is  discussed  in  depth  in 
Section  2.3.4  (Table  6). 

2.3.  Membrane  modelling 

This  section  describes  the  trans-membrane  water  flux  in  a  HT- 
PEM  FC.  The  resulting  concentration  distribution  affects  several 
sub  models  mentioned  before.  The  cross  over  of  water  is  described 
by  Fields  first  law,  which  is  derived  from  EQ:20  for  a  two  compo¬ 
nent  system.  The  two  mass  sources  represent  the  evaporation  and 
condensation  of  water  (Table  7). 


VJi  =  (po)iJ2DikVxk)  =  Ri 


(20) 


2.3.1.  Hertz— Knudsen  equation  (phase  change) 

It  is  assumed  that  the  kinetics  of  a  phase  change  can  be 
described  by  the  Hertz-Knudsen  equation  (EQ:21).  The  driving 
force  is  the  difference  in  chemical  potential  expressed  by  the 
difference  of  the  vapour  pressure  according  to  the  liquid  phase 
composition  and  the  partial  pressure  of  water  in  the  gas  phase. 
The  ‘conductivity  term’  the  square  root  expression  in  EQ:21  is 
deduced  from  the  kinetic  theory  of  gases. 

^ = ©PC(pvp|pAM  ■  ph2°)  /GtSnJ  (2i) 

The  equation  in  its  original  form  contains  two  correction 
factors,  which  in  our  case  are  simplified  to  unity  [29],  The  term 
(A/1/)pc  is  the  specific  area  of  the  vapour-liquid  boundary  within 
the  geometric  volume  of  the  catalytic  layer.  It  was  estimated  to 
be  100  [m2/m3]  and  is  assumed  to  be  constant.  pvp|PAM  is  the 
vapour  pressure  of  water  over  a  PAM  and  depends  on  tempera¬ 
ture  and  composition.  The  VLE  data  necessary  for  describing 
Pvp|PAM  were  obtained  in  our  laboratory  and  are  in  the  process  of 
publication.  pH2o  is  the  current  partial  pressure  of  water  in  the 
gas  phase. 


2.3.2.  Water  cross  over  diffusion 

During  FC  operation  water  is  produced  at  the  cathode. 
Therefore  the  environment  of  both  electrodes  in  terms  of  local 
humidity  differ  and  the  source  terms  rhpC  and  ml,5,!  will  cause  a 
concentration  gradient  within  the  membrane.  In  order  to 
describe  this  diffusive  flux  the  diffusion  coefficient  of  water 
and  phosphoric  acid  mixtures,  as  seen  in  Fig.  4,  was  obtained 
by  measurements  in  our  laboratory  and  from  the  literature 
[30,31  . 


2.3.3.  Volume  fraction  of  a  phosphoric  acid  mixture 

The  density  of  PAMs  depends  on  its  composition  and  tempera¬ 
ture.  The  mass  of  the  pure  component  phosphoric  acid  is  constant, 
therefore  a  change  of  the  density  directly  affects  the  volume 
(assumed  negligible  vapour  pressure  of  PA  up  to  200  °C  [16]).  With 
increasing  dilution  of  the  solution  its  volume  swells  and  thus  the 
volume  fraction  of  the  phosphoric  acid  mixture  within  the  catalytic 
layers  changes,  as  described  by  EQ:22  and  EQ:23. 


mPA,0  +  0-5Ampft 

Ppamwpa^cl 


(22) 


mPA,0  +  °-5AmPA 
Ppamwpa^cl 


(23) 


Am™  =  A(Ppam«pa)14<$U  (24) 

EQ.:22  shows  the  volume  fraction  of  the  phosphoric  acid  mixture 
in  the  anode.  mpA  is  the  mass  of  pure  PA,  the  index  0  refers  to  the 
moment  of  assembly.  pPAM  is  the  density  of  the  phosphoric  acid 
mixture,  wpa  is  the  mass  fraction  of  the  pure  component  PA  and  Vql 
is  the  volume  of  the  catalytic  layer.  is  the  amount  of  PA 

which  gets  pushed  from  the  membrane  into  the  electrodes  due  to 
the  swelling.  A(ppamwpa)  describes  the  change  of  the  density  of  PA 
within  the  membrane  between  operating  conditions  and  assembly. 
It  is  assumed  that  the  excess  PA  gets  distributed  to  both  electrodes 
equally,  hence  the  factor  0.5  in  EQ:22  and  23. 

Because  of  the  sum  of  all  volume  fractions  equals  unity  (EQ.:25), 
the  porosity  is  now  linked  to  the  volume  fraction  of  PAM  as  seen  in 
EQ:27.  Due  to  the  use  of  the  Bruggeman  model  the  change  of  the 
porosity  affects  different  physical  properties  summarized  in 
Table  8. 


1  =  ‘/’PAM  +  ‘/’solid  +  e 


(25) 


^  )  ^=0.5 
wpt "  Pcarb  J  VCL 


(26) 


e(T,w)  ~  1  -  </>pam(7\  wpa)  -  ‘/’solid 


(27) 


50 

Temperature  [°C] 


Fig.  4.  Diffusion  coefficient  in  the  system  water-phosphoric  acid. 
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Table  8 

Porosity  effects  in  the  CL 

Transport  of  mass  &  momentum  Porosity 

Transport  of  species  Diffusion  coefficient 

Transport  of  ions  Conductivity 

Transport  of  electrons  Conductivity 

Transport  of  heat  Conductivity 


An  important  effect  of  the  expansion  of  phosphoric  acid  is  the 
coverage  of  platinum  particles  in  the  catalytic  layer  and  thus  the 
reduction  of  ECSA,  as  explained  in  Section  2.3.4. 

2.3.4.  Electrochemical  active  surface  area 

In  Section  2.2.5  Rcat,  which  is  the  total  available  platinum  surface 
area  per  volume  electrode,  was  deduced.  However  not  all  of  that 
area  is  activated.  In  order  to  be  activated  the  platinum  particle 
needs  to  be  attached  to  an  electrically  conductive  surface  as  well  as 
to  be  in  contact  with  the  ion  conducting  phase.  As  presented  in 
Section  2.3.3  the  volume  fraction  of  the  electrolyte  depends  on 
operating  conditions.  To  account  for  the  expansion  characteristics 
of  PAM  the  relative  ECSA,  ©active,  was  introduced.  In  Fig.  5  is  shown 
how  ©active  changes  as  the  catalytic  layer  gets  filled  with  electrolyte. 
The  catalytic  area  is  considered  to  be  a  porous  medium  with  two 
distinct  pore  sizes:  One  from  the  meso  scale  and  one  from  the 
macro  scale  [32,33],  In  addition  it  is  assumed  that  PAM  is  wetting 
the  electrode  material.  When  the  electrode  gets  loaded  due  to  the 
swelling  of  the  electrolyte  the  smallest  capillaries  get  flooded  first. 
This  leads  to  a  high  penetration  of  the  electrolyte  within  the  CL  and 
wetting  of  the  surface  area  of  the  macro  pores,  as  seen  in  Fig.  6.  The 
macro  pores  covered  with  a  thin  film  of  electrolyte  are  considered 
to  be  activated.  This  condition  represents  maximum  ECSA.  With 
further  increasing  amount  of  PAM  the  activated  macro  pores  get 
filled  and  thus  inhibited  and  deactivated.  This  deactivation  results 
from  the  mass  transport  limitation,  posed  by  diffusion  through  a 
growing  PAM  layer  as  seen  in  Fig.  7.  Fig.  7  is  the  result  of  Fick’s  first 
law  of  diffusion  EQ:28  wherein  J  is  the  molar  flux,  8  is  the  diffusion 
length,  P,  is  the  partial  pressure,  H,  the  Henry  coefficient  of  the 
component  and  Qy  is  the  concentration  of  the  species  at  the  cat¬ 
alytic  surface,  which  is  assumed  to  be  zero. 

<28> 


3.  Experiment 


3.1.  Fuel  cell  test  station 

The  fuel  cell  test  station  Evaluator  C  100  by  Fuel  Con  with  the 
operating  software  FuelWork  was  used.  The  test  station  is  capable 
of  handling  loads  up  to  3  kW  and  the  heating  allows  operation  up  to 
200  °C 

3.2.  MEA  and  single  cell 

For  all  experiments  a  BASF  Celtec  P  1000  MEA  was  used.  Its 
properties  are  summarized  in  Table  9.  The  flow  fields  were  made 
out  of  carbon  with  serpentine  formed  grooves,  as  seen  in  Fig.  1. 

3.3.  Procedure 

After  a  break  in,  or  conditioning  procedure  of  100  h,  the  polar¬ 
ization  curves  were  obtained  by  using  a  voltage  guided  (potentio- 
static)  mode.  Temperature  and  pressure  sweeps  were  recorded.  The 
standard  operation  conditions  are  listed  in  Table  10. 

4.  Results  and  discussion 

4.1.  Model  validation 

In  order  to  validate  the  model  a  temperature  and  a  pressure 
sweep  was  performed.  The  results  from  the  simulations  are  shown 
together  with  corresponding  experimental  data  in  Figs.  8  and  9.  For 
the  temperature  range  the  model  predictions  are  in  good  agree¬ 
ment  with  the  experiment.  The  pressure  sweep  comparison  is  in 


For  the  validation  of  the  numerical  model  polarization  curves 
were  obtained.  Materials  and  equipment  used  are  explained  in  this 
section  whereas  the  results  are  presented  in  comparison  to  the 
numerical  studies  in  Section  4. 


Fig.  5.  Influence  of  the  PAM  expansion  on  the  ECSA:  quick  penetration  and  activation 
due  to  loaded  micro  pores  followed  by  deactivation  of  the  macro  pores. 


Fig.  7.  Growing  of  the  diffusion  barrier  with  the  film  thickness. 
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Table  9 

Technical  specifications  of  the  Celtec  P  1000  MEA  by  BASF  [37], 
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Parameter 


Electrode  material 

Cata.  material 
Pt  particle  size 
Cata.  loading  anode 
Cata.  loading  cathode 
Membrane  material  composition 
CO  tolerance 


400  pm 

4  nm  (estimated) 

1.0  mg  cm2 
PA  doped  PBI 
1000  ppm 


Table  10 

Standard  operation  conditions. 

Parameter  Value 


Anode  Cathode 


Humidity  -  18  °C  14%  r.H. 

Minimum  flow  rate  0.025  dm3  min-1  0.3dm3  mkr1 

Stoichiometric  factor  1.5  2.2 

Temperature  160  °C 


fair  agreement  as  well,  except  for  the  overpotential  region  of  the 

3  bar  polarization  curve.  All  six  U— f-curves  from  Figs.  8  and  9  were 
calculated  with  the  same  set  of  parameters,  shown  in  Table  11. 

4.2.  Sensitivity  analysis 

4.2.1.  Phosphoric  acid  content  in  the  MEA 

Fig.  10  shows  the  variation  of  the  parameter  PA  mass  in  the  MEA 
from  14.5  mg  cm-2  to  19.8  mg  cm  2.  Changing  from  14.5  to 
16.5  mg  cm~2  yields  a  slight  increase  in  performance.  Any  further 
increase  of  PA  mass  reduces  the  performance  below  the  starting 
value  and  at  19.8  mg  cm-2  the  CL  is  flooded  and  the  performance 
drops. 

4.2.2.  Humidity 

Polarization  curves  with  varying  bubbler  temperatures  (cathode 
side)  are  shown  in  Fig.  11.  With  increasing  humidity  the  perfor¬ 
mance  of  the  fuel  cell  decreases.  For  comparison  experimental  data 
is  presented  as  well.  The  experiments  were  carried  out  with  a 
different  MEA,  therefore  only  the  trends  can  be  compared. 

4.3.  PA  expansion 

The  concentration  of  PA  depends  on  the  temperature  and  on 
the  relative  humidity.  Daletou  [34]  conducted  TGAs  of  a  doped 


Current  density  [A  cm  2] 


Fig.  8.  Comparison:  simulation  and  experimental  data,  temperature  sweep. 


PBI-membrane  with  varying  partial  pressures  of  water.  A  distinct 
sorption  behaviour  of  PA  was  observed.  At  95  °C  in  dry  conditions 
the  membrane  has  the  same  weight  as  at  160  °C  with  an  enforced 
water  vapour  pressure  of  approximately  10  kPa.  Under  the 
assumption  that  only  water  evaporates  during  the  time  frame  of 
the  TGA  one  can  conclude  that  PA  has  the  same  concentration  at 
both  temperatures.  Furthermore  neglecting  the  thermal  expan¬ 
sion  of  liquids  leads  to  the  conclusion  that  the  PA  has  the  same 
volume  as  well.  This  important  characteristic  was  implemented 
into  our  model  as  described  in  Section  2.3.4.  During  fuel  cell 
operation  the  relative  humidity  depends  on  the  current  density 
(rrzH2o  —  JorrMh2o/(2F)),  therefore  on  the  operating  point.  Fig.  12 
shows  the  phosphoric  acid  redistribution  during  a  polarization 
curve.  Due  to  its  expansion  behaviour  with  increasing  dilution  PA 
moves  towards  the  catalytic  layers.  As  can  be  seen,  the  mass  of  PA 
in  the  membrane  decreases  and  the  amount  in  the  catalytic  layers 
increases.  In  sum  the  mass  of  PA  does  not  change.  In  Fig.  13  the 
temperature  dependence  of  the  expansion  behaviour  is  shown. 


Current  density  [A  cm  2] 

Fig.  10.  Sensitivity  analysis:  amount  of  phosphoric  acid  in  the  MEA. 
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Current  density  [A  cm  2] 


Fig.  11.  Polarization  curves  with  varying  humidifier  temperature  (cathode). 


Cell  voltage  [V] 

Fig.  12.  Redistribution  of  PA  within  the  catalytic  layers  during  a  simulation  run. 


The  diagram  is  divided  by  vertical  lines  to  mark  a  change  of  the 
inflow  velocity.  Due  to  minimum  flow  rates  the  inflow  velocity  for 
low  current  densities  is  not  stoichiometric.  For  the  anode  the 
minimum  flow  rate  is  reached  for  i  =  0.105  A  cm  2  and  at  the 
cathode  for  i  =  0.53  A  cm-2.  As  long  as  the  current  density  is  below 
these  values  the  inflow  velocities  are  constant  (see  Section  3.3). 
The  three  temperature  graphs,  120  °C,  160  °C,  and  200  °C,  corre¬ 
spond  to  the  temperature  sweeps  from  Section  4.1.  As  can  be  seen 
from  Fig.  13,  the  volume  fraction  of  PA  decreases  with  increasing 
temperature.  This  results  from  the  temperature  dependence  of 
the  vapour  pressure  of  water  over  PAM.  Combined  with  the  actual 
partial  pressure  of  water  the  vapour  pressure  is  the  driving  force 
of  the  evaporation,  hence  of  the  volume  change.  Therefore  the  PA 
expansion  also  depends  on  the  current  density  and  on  the  flow 
velocity,  because  both  affect  the  current  humidity. 


Current  density  [A  cm  ~2] 


Fig.  13.  Expansion  of  PAM  in  the  catalytic  layer  at  different  temperatures. 


4.4.  Flooding 

In  order  to  analyze  the  behaviour  of  a  flooded  CL  the  mass  of 
PAM  was  increased  to  20  mg  cm  2.3  Under  operation  conditions  the 
high  amount  of  PA  leads  to  a  water  uptake,  which  results  in  a 
completely  filled  CL.  The  reduced  performance  of  the  FC  is  shown  in 
Fig.  14  and  the  corresponding  volume  fraction  in  Fig.  13.  Mater  [35] 
conducted  experiments  with  excess  PAM  in  the  MEA  and  found  the 
performance  to  decrease  with  increased  PAM  content.  By 
increasing  the  amount  of  PAM  from  17.83  mg  cm-2  4  to  19.8  mg  cm2 
4  matching  results  could  be  obtained  (Fig.  15). 

4.5.  In-situ  concentration  of  phosphoric  acid 

Highly  concentrated  PA  under  STP  is  hydrophilic  and  tends  to 
absorb  water.  Because  of  the  varying  humidity  within  a  FC,  it  is 
difficult  to  estimate  the  in-situ  concentration  of  PAM.  Our  model 
predicts  the  concentrations  as  shown  in  Fig.  16.  At  0.35  V  cell 
voltage  the  mass  fraction,  wPA,  ranges  from  0.96  to  l.ll.4  With  the 
following  conceptual  model  it  is  possible  to  estimate  the  minimum 
concentration  of  PAM  in  the  membrane:  Consider  a  system  con¬ 
taining  a  gas  and  a  liquid  phase  as  seen  in  Fig.  17.  The  liquid  phase 
(2)  is  a  phosphoric  acid  water  mixture.  The  gas  phase  (1 )  is  an  inert 
gas  with  a  defined  water  partial  pressure.  The  pressure  of  the 
system  is  1  bar  and  its  temperature  is  constant  at  160  °C.  It  is 
assumed  that  at  160  °C  PA  does  not  evaporate  and  that  the  phases 
are  at  equilibrium.  Therefore  the  water  pressure  of  the  given  PAM 
equals  the  present  partial  pressure  of  water: 

Pvp|pam(T  wpa)  =  Ph2o  (29) 

The  water  vapour  pressure  of  the  a  PA-mixture  with  a  given 
temperature  determines  the  concentration  of  phosphoric  acid. 

WpA~Pvp|PAM(r)  =  Ph2o  (30) 

Assuming  the  gas  phase  consists  solely  of  water  vapour  hence  a 
partial  pressure  of  1  bar,  the  concentration  can  be  calculated  to  be 
86  wt%,  which  is  the  maximum  dilution  of  phosphoric  acid  at 
160  °C  and  1  bar. 

Compared  to  a  FC  operated  at  1  bar  and  160  °C  the  deduced 
concentration  needs  to  be  higher  because  the  partial  pressure  of 
water  is  far  less  than  the  assumed  1  bar.  Because  a  water  pressure  of 
1  bar  would  mean  the  absence  of  oxygen  and  therefore  the 
termination  of  the  FC  reaction.  In  addition,  the  membrane  in  a  FC 
faces  the  moist  cathode  as  well  as  the  much  drier  anode  gas.  The 
drier  condition  at  the  anode  side  increases  the  concentration  of  the 
PAM  further  because  of  the  evaporating  water. 

5.  Assumptions  revisited 

5.1.  Polycondensation  of  phosphoric  acid 

As  explained  in  Section  2.1  a  small  and  positive  standard 
enthalpy  of  reaction  for  the  first  polymerization  step  for  the  tem¬ 
perature  up  to  100  °C  is  known.  Because  of  the  lack  of  further  in¬ 
formations  regarding  the  temperature  dependence  of  the  chemical 
equilibrium  and  the  known  value  of  AH°(T  <  100  °C)  we  decided  to 
neglect  the  temperature  influence  and  to  use  Jameson’s  data  of  the 
composition  of  PA  for  the  discussion.  It  is  important  to  revise  that 
calculated  physical  properties  account  for  the  real  composition  due 
to  the  use  of  the  “apparent  concentration"  (see  Section  2.1). 


3  "h3po4  =  0.95@25  °C. 

4  Values  higher  than  1  are  possible  because 


1  of  polycondensation. 
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Current  density  [A  cm  2] 

Fig.  14.  Predicted  performance  of  the  flooded  FC. 


©  T  =  T" 
P'  =  P" 


Fig.  17.  Conceptual  model  for  estimating  the  minimum  concentration  of  PAM  in  a  FC. 


Table  12 

Jameson:  composition  of  the  strong  phosphoric  acids  at  15  °C  [13], 

H3P04[%]  P205  [%]  Ortho-  Pyro-  Tri-  Tetra  or  higher 

<93  67.4  100 

102  73.9  77.1  22.1  0.79 

113.7  82.4  7.3  23  19.3  50.4 


Fig.  15.  Comparison:  simulation  and  experimental  data  [Mater]  [35], 


In  Fig.  16  the  in-situ  concentration  of  PAM  is  shown  at  the 
operating  voltage  of  0.35  V.  The  maximum  concentration  is  located 
near  the  entrance  of  the  dry  hydrogen  stream.  The  mean  concen¬ 
tration  calculates  to  100  wt%.  At  this  operating  point  the  cell  pro¬ 
duces  large  amounts  of  water,  which  leads  to  a  dilution  of  PA. 
Therefore  in  Table  13  the  mean  concentration  of  PA  within  the 
membrane  for  the  whole  polarization  curve  is  summarized.  Below 
0.6  V  the  concentration  is  less  than  102  wt%  PA.  Followingjameson, 
phosphoric  acid  with  this  concentration  still  consists  of  77.1  wt% 
ortho-PA  (Table  12).  Especially  at  lower  operating  voltages  and 
higher  temperatures  polymerization  might  occur.  It  seems  neces¬ 
sary  to  further  investigate  the  effects  of  polycondensation 
regarding  the  physical  properties  of  highly  concentrated  phos¬ 
phoric  acid  containing  poly-phosphoric  acids. 


5.2.  Volume  swelling  of  the  PBI-membrane 
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The  thickness  of  the  membrane  depends  on  the  acid  content 
[17],  Fig.  13  shows  the  amount  of  phosphoric  acid  in  the  membrane 
during  a  polarization  curve.  As  described  in  Table  13,  20%  of  the 
mass  of  phosphoric  acid  leave  the  membrane  towards  the  CL.  It  is  to 
be  expected  that  this  amount  will  change  the  thickness  of  the 
membrane,  which  until  now  is  not  accounted  for. 

6.  Conclusion 

A  2D  FEM  model  of  a  HT-PEM  fuel  cell  was  built  and  enhanced 
by  adding  the  vapour  liquid  equilibrium  and  evaporation  kinetics 
into  the  coupled  system.  Thereby  the  substance  distribution  of 
water  respectively  phosphoric  acid  was  computed  and  used  for 


Table  13 

Cell  voltage,  current  density,  mass  density  of  PA  in  the  membrane  and  its  mean 
concentration  during  a  polarization  curve  simulation. 


ifCeii  [V]  MAcmj  mPA,Mem  [mg  cm-2]  0  6,PA[wt%] 


0.9 

0.77 

0.72 

0.70 


0.00 

0.01 

0.03 

0.07 

0.10 

0.21 

0.34 

0.49 

0.67 

0.88 

1.10 


11.41 

11.35 

11.32 

11.30 


114 


Fig.  16.  In-situ  concentration  PAM  0.35  V. 
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determination  of  physical  properties  while  solving  the  model.  The 
in  situ  concentration  of  phosphoric  acid  was  found  to  be  in  the 
range  of  96  to  111  wt%. 

Polarization  curves  for  different  temperature  and  pressure 
sweeps  was  simulated  and  compared  to  the  experimental  data.  The 
model  showed  consistency  and  was  used  to  investigate  the  sorption 
behaviour  further.  The  simulation  is  able  to  predict  the  effects  of 
excess  acid  as  investigated  by  Liu  [35],  The  following  conclusion  can 
be  drawn: 

Humidity  strongly  affects  the  performance  of  high  temperature 
fuel  cells.  The  growing  of  the  acid  film  layer  due  to  sorption  can 
cause  severe  mass  transport  limitations.  This  effect  can  be  reduced 
by  increasing  the  temperature,  which  leads  to  an  increased  vapour 
pressure  of  water  and  a  higher  concentration.  Reducing  the  partial 
pressure  of  water  in  the  gas  phase  by  lowering  the  current  density 
resp.  the  water  production  has  the  same  effect.  As  the  intended 
application  for  high  temperature  fuel  cells  is  the  operation  with  an 
upstream  steam  reformer  and  therefore  high  humidity  conditions, 
the  investigation  of  the  water  balance  is  of  great  importance. 

Regarding  the  effects  of  humidification  and  flooding  it  is 
important  to  point  out  that  the  condition  of  the  MEA  is  crucial  for 
the  result.  Obviously  when  diluted  the  ionic  conductivity  of  PA 
increases.5  Many  publications  observed  the  effects  of  humidifica¬ 
tion  with  differing  system  responses  [34-36],  The  discrepancy 
whether  increasing  the  humidity  of  a  HT  fuel  cell  increases  the 
performance  or  not  depends  on  the  content  of  PA  in  the  catalytic 
layer.  An  increase  in  performance  as  observed  by  Daletou  [34] 
might,  in  addition  to  the  increased  ionic  conductivity,  be  the  ef¬ 
fect  of  a  lack  of  PA  within  the  CL.  With  increasing  expansion  of  PA 
due  to  sorption  of  water  more  catalytic  particle  are  connected  to  the 
electrolyte  phase  and  the  ECSA  grows. 

The  opposite  effect  as  described  by  Liu  [35]  and  Lobato  [36]  and 
also  observed  by  our  bachelor  student  M.  Bahr  (Fig.  11 )  might  be 
caused  by  a  higher  amount  of  PA  in  the  CL  which  leads  to  decent 
performance  under  unhumidified  conditions.  But  the  expanded  PA 
during  humidified  operation  will  pose  a  deactivation  of  ECSA  and 
hence  a  decrease  in  performance. 

This  work  is  intended  to  be  the  basis  for  future  publications 
regarding  degradation  due  to  water  stress  and  3D  simulations  for 
determining  the  water  discharge  rate  at  the  anode.  Furthermore  we 
are  preparing  the  publication  of  our  VLE  measurements  regarding 
the  vapour  pressure  of  phosphoric  acid  and  the  vapour  pressure  of 
water  in  the  presence  of  phosphoric  acid. 
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